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1. INTRODUCTION 



Since the discovery of "Dirac equation" [Ref.l], there have been many years of 
development of relativistic quantum theories on Dirac fields. It is known that 
the Dirac equation is "metrical" , in which spinor has more than one component. 
There are many references on spinors, like the book by Dirac [Rcf.2] and the 
standard textbook by Bjorken and Drell [Ref.3]. The existence of spin in the 
Dirac equation, is a natural outcome of the unification of special relativity and 
quantum mechanics, and is one of the greatest triumphs of the Dirac theory. 
The relationship between spinor and spacetime has long been an interesting 
research subject, and there are many different approaches with a huge amount 
of publications along this line. But I have always been thinking about tackling 
these issues in a more straightforward and easy-to-follow kind of fashion. This 
paper is written exactly in this way. Hope it can bring some new results and 
insights to the development of Dirac quantum field theory. 

Given those well-known results on spinors, I intend to explore the possibility 
of dealing with spin and spacetime in a simply unified fashion. The main idea 
is make use of a minimal complete set of variables with both spin and space- 
time degrees of freedom, which I shall call "common variable" (see Jin [Ref.4]). 
It is then straightforward to construct a class of "common operators" such as 
"common momentum" out of this variable. I shall give this space, established 
directly on the idea of "common variable" , the name "common space" for clarity 
There are physical as well as mathematical reasons to construct such a "common 
space" . Physically, spin, an "intrinsic" property of the electron or other spin- 
half fcrmions, appears naturally with respect to the "external" Lorentz trans- 
formations in Minkowski spacetime. There is always a sense of correspondence 
between spin and spacetime. Mathematically, this "common space" possesses a 
complete set of degrees of freedom in both spin and spacetime, and transforms 
in an "invariant" way, with spin and spacetime each transforming in a "covari- 
ant" way with respect to the Lorentz transformations. It is a simple union of 
spin and spacetime, in the sense of special relativity. 

One of the most surprising results I come up with, is that the "common 
momentum" and "common variable" of a massive fermion, constructed in this 
"common space" , are exactly the kind of linear self-adjoint "doubly strict plus- 
operators" in an indefinite inner product space, as sought-after and discussed 
in detail in an excellent mathematical book by Bognar [Rcf.5]. There is a long 
history of research in indefinite inner product spaces, starting from Dirac [Ref.6], 
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then Pauli [Re. 7], Heisenberg [Refs.8,9] and many other researchers. The idea 
is to find a finite representation like the Dirac spinor representation for the 
homogeneous Lorentz group including non-compact Lorentz boosts, which can 
only be given in indefinite inner product space. There are some unresolved 
issues in this field, such as the definitions of positive probability, vacuum state, 
etc. I will provide my answers to these questions in this paper. I will also prove 
an uncertainty relation between fcrmion mass and proper time, by utilizing 
the commutation relation of "common momentum" and "common variable". 
Finally, I will present Dirac quantum field theory in the notion of "common 
space" with the introduction of so-called Lagrangian operators, which need to be 
self-adjoint in the indefinite inner product space to preserve Lorentz invariance. 
Nonlinear QED equations and Lagrangians [Ref.4] will be given as an example 
to demonstrate this approach. 

2. COMMON SPACE 

To begin with, we define a common variable in terms of Dirac matrices 7^ and 
spacetime coordinates: 

n = -fX^ (2.1) 

where 7^ are expressed in the standard representation [Ref.3]. The derivative 
with respect to this variable can be deduced as 

^ = (^q)^ = ^- ( 2 - 2 ) 
Hence a common momentum can be introduced as 

P u = id u = -fP ll . (2.3) 

The natural units are to be used throughout, unless specified otherwise. It is 
straightforward to show the following commutation relation 

[Sl,P u ] = -U, (2.4) 

by the anticommutation relations of 7 M : {7^, 7"} = 2g flu , and the commutation 
relations of X^ and P v : [X^, P v ] = —ig^, where g^ u or g^ v are spacetime metric 
elements with Minkowski signature (1, —1, —1, —1). Here we have utilized the 
commutation relations of time and energy as well as space and momentum, 
discussed by Stuecklbcrg [Rcf.10], Horwitz and Piron [Refill]. 
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In common space, the Lorentz transformation L u is a direct product of a 
spinor one L s and a coordinate one L c : L u = L S L C . Since momentum P M are 
covariant in spacetime 

LcP^L- 1 = a^P v , (2.5) 
and Dirac matrices 7 M are contravariant in spinor space 

L^L- 1 = aW, (2.6) 

it is easy to check that 

LujPujL^ 1 = Pv, (2-7) 

namely, common momentum is Lorentz invariant, and so are any other common 
operators to be introduced in common space. 
Now for any two common operators 

Y u = y%, (2.8a) 

Z u = 7^, (2.86) 

we define a scalar product that eliminates spin components and provides a scalar 
in spacetime only: 

Y u -Z u = Y^. (2.9) 

A simple example is P 2 — P^ ■ P u = P^P^ . 

Since the common momentum squared of a free particle with mass to is 
constant P 2 = E 2 — p 2 = to 2 , we have a wave equation of the second-order 
Klein- Gordon type 

P 2 *(ft) = mH(n), (2.10) 
and two wave equations of the first-order Dirac type 

P W V{Q) = ±m*(fi). (2.11) 

In this simple case, P w is a physical observable with eigenvalues ±to. 

The common momentum and common variable so-defined are not hermitian 
but instead pseudohermitian: 

7 o P( t 7 o = p ^ (2 12fl) 

7°fiV = n (2.12b) 

Consequently they are not well-defined in positive-definite Hilbcrt space where 
all self-adjoint operators are hermitian. To make them physical observables, we 
need to find a better way. 
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3. INDEFINITE INNER PRODUCT SPACE 



Back in 1942, Dirac [Ref.6] initiated the idea of implementing an indefinite inner 
product space. This idea was subsequently explored by Pauli [Ref.7] who found 
a self-adjointness requirement: the operators of physical observables in such a 
space be pseudohermitian rather than hermitian. Hciscnbcrg [Refs.8,9] later 
tried to implement this idea in his unified field theory as well. 

In a positive-definite Hilbcrt space, the invariance group is restricted to the 
compact unitary group including rotations and inversions, while the Lorentz 
boosts which are non-compact and non-unitary (though unimodular) arc not 
included. This situation can be clearly stated according to Heisenberg [Ref.8]: 
"a finite representation of a non-compact group can be given only in space with 
indefinite metric" . 

Since the early work by these pioneers, there has been a lot of progress in 
this field. There is an excellent book by Bognar [Rcf.5] on the mathematical as- 
pects of indefinite inner product spaces. In this book, emphasis has been put on 
Krein spaces, the most important type of inner product spaces, which stand for 
non-degenerate, decomposable, and complete spaces. Interested readers should 
study it for detailed mathematical treatments. We shall only utilize some im- 
portant definitions and results from this book, which are directly relevent to 
our approach. In the rest of this paper, we shall assume all our inner product 
spaces are Krein spaces. 

For our purpose, we define an indefinite inner product space (HPS) for Dirac 
fields, which is an indefinite Krein space, with the following inner product 

< ip,<p >= J drtjxj), (3.1) 

where dr = d x, and ip = ip'-f° is the adjoint of ip. It can be divided into three 
subspaces with zero, positive and negative norms respectively 

Z = {ip\ < 1>,1> > = 0}, (3.2a) 

V= {V| < ip,ip > > 0}, (3.26) 

Af= {ip\ < ipA > < 0}. (3.2c) 

A simple example for a zero norm subspace is ip = cl G Z, a four-component 
column unit vector / multiplied by a constant c. It is the negative norm subspace 
that causes trouble and confusion. 
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In such an HPS, the common momentum P u and common variable f2 are 
well-defined physical observablcs since they are both self-adjoint operators, 
namely, their expectation values are real 

<P^>- [ ^ ± -<P^>\ 0.3a) 
J dripip 

<n>= =< n >t . (3.36) 

J dripip 

To make the above integrals converge, we need to assume certain boundary 
conditions and let ip £ Z. 

Mathematically, a symmetric operator in HPS obeys: 

< V,0*</> >=< Oip,<p >=< ip,0(f> >, (3.4) 

namely its adjoint is 

O* =7 °Ot 7 = o (3.5) 

A self-adjoint operator in HPS is a symmetric operator with dense domain. We 
shall assume all our symmetric operators satisfying (3.5) have dense domain 
in HPS, and simply call them self-adjoint operators. In this sense, common 
momentum and common variable are self-adjoint by (2.12a) and (2.12b). Note 
that if O commutes with 7 : [0, 7 ] = 0, then it is hermitian. 

A transformation L keeping the indefinite inner product invariant yields: 

< Lip, L<j> >=< V, <t> > ■ (3-6) 

This is true for almost all the homogeneous Lorentz transformations including 
non-compact Lorentz boosts [Ref.3], which satisfy 

^Ltj°L u = I. (3.7) 

Note that if a Lorentz transformation commutes with 7 : [L w ,7°] = 0, then it 
is a unitary transformation such as a rotation or a space inversion. It is special, 
however, for unitary time inversion [Ref.4] that satisfies 

7 Tt 7 T u , = -/. (3.8) 

Given that the time integral by dt reverses sign under T u , unitary time in- 
version also renders indefinite inner product invariant in (3.6). These results 
show that the so-defined indefinite inner product space provides a finite spinor 
representaion for the whole non-compact homogeneous Lorentz group. 
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A major concern about the legitimacy of HPS is whether one can come up 
with a valid interpretation of positive probability [Ref.9]. If simply defining a 
probability density as 

po(x) =#r)V>(z), (3.9) 

one would get into trouble for po(x) < when ip(x) <E AT. To avoid this difficulty, 
we need to redefine the probability density in such a way that total probability 
adds up to one: 

I**- 1 - < 3 - io) 

With this definition, we now have a normalized positive probability density 
p(x) > no matter ijj(x) e V or M unless ip(x) E Z. This new definition is 
more natural, in the sense that it provides not only a generalization from non- 
relativistic to relativistic quantum mechanics by preserving Lorentz invariance, 
but also a reasonable physical interpretation in a statistical fashion. 

In Chapter II of Ref.5, there are some detailed discussions about so-called 
fundamental symmetry J, which are quite analogous to the way we define our 
normalized positive probability. One can set J = P + — P~ with the fundamental 
projectors P + and P~ onto the positive and negative subspaces respectively. 
Then one can define a J-inner product 

(4>,x)j = (W,x) = (V> + ,x + ) - (V",X"), (3-11) 

which is positive-definite, namely, (ip,ip)j > 0. 

4. DOUBLY STRICT PLUS-OPERATOR 

We need some definitions out of Bognar's book, for subsequent discussions. A 
linear operator in HPS is called a "plus-operator" , if it is defined everywhere in 
HPS and carries non-negative vectors into non-negative ones. A "plus-operator" 
is called a "strict plus-operator" , if it carries positive vectors into positive ones. 
A "strict plus-operator" is called a "doubly strict plus-operator" , if its adjoint 
is also a "strict plus-operator" . 

For a linear operator T in HPS, define 

fi{T) = infimum <lpMl>=1 < Tip, Tip > ■ (4.1) 

Then it turns out that T is a plus-operator if //(T) > 0; a strict plus-operator if 
/z(T) > 0; and a doubly strict plus-operator if (a(T) > and fi(T*) > where 
T* is the adjoint of T. 
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For a massive fermion, it happens that its common momentum and common 
variable are both "time-like" , namely 

< H >= <P f >™> > 0, (4.2a) 

~ < nip, nip > 

< n 2 >= T—r- — > °> ( 4 - 2& ) 

<ip,ip> 

therefore they are "strict plus-operators". Furthermore, they are self-adjoint, as 
mentioned in the last section, so their adjoints are also "strict plus-operators" . 
We thus come to a very important conclusion: the common momentum and 
common variable of this massive fermion fall into the category of "doubly strict 
plus-operators" . 

A nice feature about this kind of operators, denoted by O now, is that the 
three subspaces are separable 

0{Z;V;N) e Z;V;M, (4.3) 

and it is easy to see from (4.2a) and (4.2b) that < P u ip, P u ip > and < flip, Clip > 
have the same sign as < ip, ip >. This can be used to prove a more general result 
(see Lemma II. 8. 7 of Ref.5): 

< Otp, OtP' >= n{0) < ip, ip' > {ip, ip' € HPS) (4.4) 

where fi{0) > for doubly strict plus-operator O. 

If < ipo,ipo >= in the zero norm subspace, then < Otpo,Oipo >— by 
(4.4). Let ipi = Otpo, then < tpi,ipi >= and < tp + rpi , tpo + ipi >= 
lead to Re < ipo,ipi >— 0. Here < ip$,OipQ > cannot be complex for O 
is self-adjoint. It has to be zero. In any case, < tp ,OtpQ >= whenever 
< tp , ip a >— 0. As a result, no physical observables can be measured in such a 
zero norm state. Heisenberg once called it "ghost" state [Ref.8]. Here we prefer 
to call it "vacuum" state. 

The conventional vacuum state is defined by 

PjXo>=0. (4.5) 

Let ip = xo ® I, we get 

Po#o>=0, (4.6) 

and < P^ipo, Puipo >— 0. From (4.4) we have < V'OjV'o >= 0. Hence all ipo 
which have zero norm become our new vacuum state in HPS, related to the 
conventional vacuum states xo in the above sense. 
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5. UNCERTAINTY RELATION 



The simplest type of plus-operators can be expressed by a doubly strict plus- 
operator O multiplied by a constant number c, since it is obvious (for tp Z) 

<cOtp,cOtp> c*c<Oip,Oip> „ 

; — ; = ; — ; ^ °- I 5 - 1 ) 

< v, i> > <ip,ip> 

The operators cP u and cfi are both plus-operators though they may not be self- 
adjoint or doubly strict when constant c is complex or simply zero. Furthermore 
any operator of the type c\P^ + C2$l is a plus-operator by the superposition of 
indefinite inner product subspaces. This can be stated as follows: if tp G V, then 
CiPu,tl> (and c 2 0^) G V+Z, we have (c\P w + C2fl)ip G V+Z; similarly if tp G J\f, 
we have (ciP w + C2^)ip G Af+ Z. In general, so long as ip £ Z, we always have 
an inequality for any constant numbers c\ and c 2 : 

< (cifi, + c 2 fl)tp, (ciP u + c 2 n)ip > >Q 
<ip,ip> ~ 

Let us define a norm ratio as a function of a real parameter £ 

< (P, - igj^ (P, - *w > 

= < P 2 > +i£ < [fi, P u ] > +£ 2 < f> 2 > . (5.3) 

To have nonzero physical observables in any non-vacuum state, from (5.2) we 
know > for any real parameter £. That is to say 

< i[fi, P,,,] > 2 -4 < P 2 >< ft 2 >< 0. (5.4) 

The same argument can be applied to the fluctuations of common momentum 
and common variable: 

SP U = P u - < P„ >, (5.5a) 

SQ = fl— < O > . (5.56) 

Noting that from (2.4) 

< [SQ, SP^} >=< [fi, P u ] >= -4i, (5.6) 
we have the following uncertainty inequality 

AP^Afi > 2ft, (5.7) 
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where 

AP W = y/< {5P„Y >, (5.8a) 

Afl = V< (^) 2 >• (5.86) 

For a free particle AP W = Am and Af2 = At, we have an uncertainty relation 
between its mass and proper time 

AmAr > 2h. (5.9) 

So there exists uncertainty in measuring particle mass, given the uncertainty 
in locating its spacetime position. This uncertainty relation was also discussed 
early by Arshansky and Horwitz [Ref.12]. 

There is an alternative way to prove inequality (5.4) by utilizing Schwarz 
inequality in a semi-definite inner product space (see Lemma 1.2.2 of Ref.5): 

(x,x)(y,y)>\(x,y)\ 2 . (5.10) 

Now we define a positive-definite inner product 

(Qip,Qip) = , (5.11) 

< ip, ip > 

for any plus-operators Q and Q' with ip £ Z. We know and «f2 are both 
strict plus-operators, namely, 

(P^,P^) >0, (5.12a) 

(iflip,iQtp)>0. (5.12b) 
From the Schwarz inequality (5.10) we get 

(p^-Pu,V0W,^V>) > \(P u il>,ini>)\ 2 

> \Re{P^,i^)\ 2 = i|(?M[ft,P,]V)| 2 . (5.13) 
This is the same as (5.4). It is similar to prove (5.7) by the Schwarz inequality. 

6. QUANTUM FIELD THEORY 

In quantum field theory, one normally starts from a Lagrangian, then applies 
variational principles to derive field equations. Using our notations, we write 
an action as an integral of a Lagrangian density C 

W= f duC(il>,d u il>), (6.1) 
Jr 
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where du> = \J (dVL) 2 = d A x is a measure in spacetime, and R is a four- 
dimensional region with a three-dimensional boundary B that can be either 
fixed or not fixed (see Barut [Ref.13]). Here we discuss fixed boundary only. 

To render the whole approach manifestly Lorentz invariant, we introduce a 
Lagrangian operator C in the indefinite inner product space, so that the action 
can be written as 



W= / duj^{x)Ctl){x) =< Lip > . (6.2) 
Jr 

Since the action is real- valued, the Lagrangian operator ought to be self-adjoint, 
namely, 

C = 7°£%° = L (6.3) 

This is a stringent condition in our approach, in contrast to the conventional 
approach of choosing Lagrangians with much arbitrariness. 

The variation of the action (6.1) with respect to field variables is given by 

dC 

dduitp 



f dC 

sw = I^W*" 



In the first term, the common operator indicates that the derivatives are 
acting on the left while the order of matrix multiplication is not altered. In the 
second term, we define 

Its integral can be changed to a surface integral by Gauss' theorem 

dC f 

du[d UJ {——-5'4)))= / dcTu-Foj, (6.6) 



where the surface common operator is defined by 

^ = 7*V ( 6 - 7 ) 
with (T M being a surface four-vector, and 

dC 

= aw 6 *' (6 ' 8) 

Considering matrix multiplication, we rewrite the right-hand side of (6.6) 

dv u ■ T u = J (^L.)d* u (<ty). (6.9) 
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Then let SW = and field vanish on the surface, we get the following 
equation of motion 

§-O a "= a (6 - 10) 

The same argument is true when ip is the field variable, which leads to 

-*,(-*£_) = (). (6.11) 

dip d(ipd^) 

For a simple free particle, we write a Lagrangian density 

£ = ipid^ip — mipip . (6-12) 

From (6.11), we get the free Dirac equation 

id w ip - mip = 0. (6.13) 

In general, for any c-number potential V(ip, ip) not involving the derivatives of 
field variables, we introduce a self-adjoint Lagrangian operator satisfying (6.3) 

£ = id w -m-V(ip,ip), (6.14) 

and write a Lagrangian density 

£ = ipid^ip — mipip — ipV(ip, ip)ip. (6.15) 

From (6.11), the equation of motion turns out to be 

—dV 

id^ip — mip — Vip — ip^=ip = 0. (6.16) 

dip 

In the case of electrodynamical interaction, we identify 

F(V^) = ^7^, (6.17) 
to obtain a nonlinear equation of motion [Ref.4]: 

id w ip - mip - e J^A^ip = 0. (6.18) 
Here four-current J M = ip^/^ip satisfies a continuity equation 

8^ = 0, (6.19) 
and four-potential is generated by an external four-current 

DA^ = qj; xt , (6.20) 
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under Lorentz gauge condition 

d^Ap = 0. (6.21) 

Equations (6.18)-(6.21) are what we need to establish a nonlinear QED in which 
Lorentz gauge is unique and inherent [Ref.4]. 

Defining a common current and a common potential 

Jo, = (6.22a) 

A, - (6.226) 
we rewrite the above set of equations by utilizing scalar product (2.9) 

iduil> - mip - e(Juj ■ A w )ip = 0, (6.23) 

B u ■ J u = 0, (6.24) 
DA, = q.PJ\ (6.25) 
d u ■ A u = 0, (6.26) 

where □ — ■ — d^d^. 

This theory is also invariant under the translations in the inhomogeneous 
Lorentz group (or Poincare group), with common momentum as the generator. 
Given a state vector ip(^o) a t a particular common variable Oo, similar to the 
Schrodinger picture 

id^(n ) = p^(n Q ), (6.27) 

we have a state vector at another common variable ft + £1 

iP(n + fi„) = [exp(-iP u • n)]V(fi ). (6.28) 

We may also have a dchnition in analogy to the Heisenberg picture. Given a 
physical common operator O u (£lo) at Oo, we define this operator at another 
common variable il + SIq by 

o u (n + n ) = [ex P (iP w • si)]o u (iio)[exp(-iP u • n)]. (6.29) 

Then we have an equation of motion 

d u -O u =i[P u ,O u ]. (6.30) 

This implies that common operator O u is conserved with respect to an evolution 
of f2, if it commutes with common momentum P^. 
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In the case of electrodynamical interaction, comparing (6.23) with (6.27) we 
obtain its common momentum 

P w = m + eJ w -A w , (6.31) 

where mass and charge are physical observables as part of common momentum. 
This provides an alternative way to unravel the mystery of mass and charge. It 
is easy to verify 

[P u ,J u ]=0, (6.32) 
to obtain conservation equation (6.24), namely continuity equation (6.19). 
We now define a conjugate common momentum of field ip 

dC 

n=— — . (6.33) 

By the action principle, a generator that provides conservation laws in integral 
form, can be defined as follows 

/ da UJ -T w = J n(da u )5il>, (6.34) 

with a being a three-dimensional surface. For electrodynamical interaction 
(6.17), Lagrangian (6.15) gives a conjugate common momentum 

n = kj, (6.35) 

and a generator 

/ da u -J r u = / da^iipi^Sip), (6.36) 

J (7 J (7 

leading to the continuity equation (6.19) in differential form as well. 

For any arbitrary separation, the anticommutation relations of Dirac field 
operators can be written as 

jr/>(fi),V(fi')} = °. (6.37a) 
{^(fi),^(fi')}=0, (6.376) 
j^ffi),^')} = -iS(n - n'), (6.37c) 

with 

S(fl) = f d Pu —^— exp(»P w • O), (6.38) 

(27r) 4 J Pu+m 

where dp^ = y/ (dP u ) 2 = d 4 p is a measure in energy-momentum space. The 

above anticommutation relations of canonical quantization always hold true, so 

long as the interaction potential in the Lagrangian (6.15) does not involve the 

derivatives of field variables. Note that there are no derivatives in our nonlinear 

QED interaction potential (6.17). 
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7. REMARKS 



This paper presents a new way to develop relativistic quantum theory on Dirac 
fields. The emphasis has been put on the construction of a "common space" 
which is actually a simple union of "intrinsic" spin and "external" spacetimc. 
A "common variable" in common space is a minimal complete set of variables 
in spin and spacetimc. A "common momentum" in common space is defined in 
terms of the derivative with respect to "common variable" . Along this line, a 
whole class of "common operators" can be defined in terms of Dirac matrices and 
spacetime operators, which are manifestly Lorentz invariant. These common 
operators are not hermitian but instead pseudohermitian, and are well-defined 
physical obscrvables in an indefinite inner product space. 

In nonrclativistic quantum mechanics, the positive-definite Hilbert space has 
been widely used in Euclidean three-space with time being a parameter. While 
in relativistic quantum mechanics, spacetime is Minkowskian with time being 
an independent coordinate. The invariance group of 4-dimensional Minkowski 
spacetime is the non-compact Lorentz group including Lorentz boosts. For 
such a non-compact invariance group, the finite representation like the Dirac 
spinor representation discussed exclusively here in this paper, can only be given 
in indefinite inner product space. To resolve one of the major contradictions 
in this approach, we have defined a normalized positive probability, which is 
Lorentz invariant as well. 

This paper furthur demonstrates that common momentum and common 
variable of a massive fermion are linear self-adjoint "doubly strict plus-operators" 
in the indefinite inner product space defined for Dirac fields. This critical result 
leads to the introduction of "vacuum state" in the zero norm subspace; to the 
uncertainty relation between fermion mass and proper time; and to the better 
understanding of Lorentz invariant physical quantities such as mass and charge. 
Recasting Dirac quantum field theory in the notion of common space, we have 
introduced self-adjoint Lagrangian operators in order to write the full action 
integral in the indefinite inner product space. By variational principles, we can 
then deduce quantum fermion field equations for different interactions, such as 
nonlinear quantum electrodynamical interaction. 

We may just as well generalize what we have done here into curved spacetime. 
First we may define an indefinite inner product space for Dirac fields with a 
generalized measure. Then we may define so-called common operators of spin 
and spacetime, such as common momentum by utilizing spin affine connections, 
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and show that these common operators are symmetric or self-adjoint in the 
indefinite inner product space. In any static spacetime with time-independent 
metric, which includes Minkowski spacetime as a special case, we can utilize 
on-mass-shell projections to project fermion Fock space onto the positive and 
negative energy subspaces, and quantize free Dirac field operators in terms of 
the creation and annihilation operators in both positive and negative energy 
subspaces, as demonstrated in another paper by Jin [Ref.14]. 

For interacting Dirac fields, quantization procedures become too complex, 
if not impossible. Once more complicated interactions are involved, nonlinear 
effects tend to be dominant, and particle creation and annihilation are not really 
clear-cut during the processes. Many unstable intermediate particles appear in 
high energy physics experiments. Certain exotic states of matter may form at 
short length scale as discussed by CP. Kouropoulos [Ref.15]. We eventually 
have to go back to some nonlinear field equations and try to find better ways to 
solve them in a nonperturbative but physically and mathematically attainable 
fashion. This is the kind of challenge we are still facing in solving nonlinear 
quantum field theories like nonlinear QED [Ref.4]. 
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